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1 Introduction 



Recently, the so called PT- symmetric quantum mechanics has attracted wide attention [1]. In Bender 
and Boettcher's work in 1998 [1], it was shown that the class of non-Hermitian, Hamilton operators 
such as H = p 2 + x 2 (ix) e (e > 0) has a real spectrum due to its PT-symmetry where P and T are 
the parity and time reversal operators respectively [2] . Exact solution of the Schrodinger equation 
for various potentials which are complex are generally of interest. It is also known that PT-symmetry 
does not necessarily lead to completely real spectrum, and an extensive kind of potentials of real 
or complex form are being faced with in various fields of physics. In particular, the spectrum of 
the Hamiltonian is real if PT-symmetry is not spontaneously broken. Recently, Mostafazadeh has 
generalized PT symmetry by pseudo-Hermiticity [3]. In fact, a Hamiltonian of this type is said to 
be rj- pseudo Hermitian if H + = rjHrj -1 , where + denotes the operator of adjoint. In [4] new class 
of non-Hermitian Hamiltonians with real spectra was proposed which are obtained using pseudo- 
symmetry. Moreover, completeness and orthonormality conditions for eigenstates of such potentials 
are proposed [5]. In the study of PT-invariant potentials various techniques from a great variety of 
quantum mechanical fields have been applied such as variational methods, numerical approaches, 
Fourier analysis, semi-classical estimates, quantum field theory and Lie group theoretical approaches 
[5-14]. In addition, PT-symmetric and non-PT symmetric and also non-Hermitian potential cases 
such as oscillator type potentials [15], a variety of potentials within the framework of SUSYQM 
[16-19], exponential type screened potentials [20], quasi/conditionally exactly solvable ones [21], 
PT-symmetric and non-PT symmetric and also non-Hermitian potential cases within the framework 
of SUSYQM via Hamiltonian Hierarchy Method [22] and some others are studied [23-25]. 
The QHJ formalism, which is a formulation of quantum mechanics was investigated as a theory re- 
lated to the classical transformation theory [26-27]. It was formulated by Leacock and Padgett in 
1983 [28-29]. Within the Quantum Hamiltonian Jacobi approach (QHJ), which follows classical me- 
chanics, not only the the energy spectrum of exactly solvable (ES) and quasi-exactly solvable (QES) 
models in quantum mechanics but eigenfunctions can also be determined [30-36]. The advantage 
of this method is that it is possible to determine the energy eigenvalues without having to solve 
for the eigenfunctions. In this formalism, singularity structure of the quantum momentum function 
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p(x) which is a quantum analog of classical momentum function p c determines the eigenvalues of 
the Hamiltonian. An exact quantization condition is formulated as a contour integral, representing 
the quantum action variable, in the complex plane. The quantization condition leads to the number 
nodes of the wave function. The wavefunction is related to the quantum momentum function (QMF). 
The equation satisfied by the QMF is a non-linear differential equation, called as quantum Hamilton- 
Jacobi equation. There is a boundary condition in the limit QMF which is used to determine physically 
acceptable solutions for the QMF [29-37]. In the applications, Ranjani and her collaborators applied 
the QHJ formalism, to Hamiltonians with Khare-Mandal potential and Scarf potential, characterized 
by discrete parity and time reversal (PT) symmetries [31]. 

The purpose of the present work has been to apply the QHJ formalism to the Hamiltonian in one 
dimension with non-hermitian exponential type potentials in order to see possible singularities for the 
QMF which determine eigenvalues and convenient eigenf unctions. 

The organization of the paper is as follows. In Sec. U, we briefly introduce the Quantum Hamilton- 
Jacobi formalism. In Sec. Ill and Resolutions of PT-/non-PT-symmetric and non-Hermitian forms of 
the well-known potentials are presented by using QHJ method. We discuss the results in Sec. V. 

2 Quantum Hamilton- Jacobi Formalism 

In quantum theory, one assumes that function W(x, E) satisfies (2m = 1) [33], 

~dW(x,EY" 



. h d 2 W(x,E) 
dx 2 



dx 

Eq. (1) will be called as the QHJ equation. The momentum function 



(E-V(x)) (1) 



, v dW(x,E) 

will be called as the QMF. In the limit U — > 0, the QHJ equation turns into the classical Hamilton- 
Jacobi equation. Then, QMF turns into the classical momentum function in the % — > limit: 



p(x, E) - p c (x, E) = ^E-V{x) 
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(3) 



In terms of p(x, E) the QHJ equation, Eq.(l) can be written as 



p 2 {x, E) - ihp'{x, E)-[E- V(x)} = 0. (4) 

Leacock and Padgett [28,29] proposed using the following quantization condition for the bound states 
in order to obtain eigenvalues. C is a contour that encloses the moving poles between the classical 
turning points and the integral 

J(E) = J- / p(x) dx (5) 

271 JC 

is called the quantum action variable. More details can be found in the paper of Bhalla et all [30-37]. 
Then 



J = nh = J(E) (6) 

gives the exact energy eigenvalues (n — 0, 1, 2, ...) [28-33]. Leacock [28,29] defines the wave function 
in order to connect QHJ equation to the Schodinger equation, 



ip(x, E) = exp 



l -W{x,E) 



(7) 



hence ip(x, E) satisfies the Schodinger equation and the physical boundary conditions. The quantiza- 
tion condition becomes [33] 

j> p(x, E)dx = 2ni ^2(Res) k = nh (8) 
Jc k 

where Y,k(R es )k is sum of the residues. In the QHJ equation, if V(x) has a singular point, p(x, E) 
will also have singular point in that zone [33]. These singularities are known as fixed singular points 
which are energy independent. Other types of singular points are the moving singular points. They 
can only be poles with residue —%h. Suppose 6^0 then, moving singularities are in the form of 
[30-37] 



p{x,E) 



+ ... 



(9) 



(x - X ) 

in the QHJ equation. If the potential is not singular at x = x then r must be equal to one and b = —ih 
[33]. 



3 Generalized Morse Potential 



The generalized Morse potential is given by [19] 



V(x) = V ie - 2ax - V 2 e~ ax 



(10) 



In order to apply QHJ method, we write the potential relation in Eq.(4) (% = 2m = 1) 



p 2 - ip -[E- V ie ~ 2ax + V 2 e~ ax } = 



(ID 



Substitution of the transformation of y = \fV[e ax in Eq.(l 1) gives: 



p 2 (y,E) +iayp'(y,E) 



2 ^2 

E - y 2 + -^=y 



V 1 ' 



= 



(12) 



Define p = iay<f> and x = + ^ in order to transform Eq.(l 1) into a Riccati type differential equation 



as, 



x + r + -r^ + 



1 



Ay 2 a 2 y 2 



E-y 2 + ^y 



(13) 



As it can be seen from Eq. (13), x has a pole only at y = and for y = define x as, 



X = h a + a x y 

y 

Substitute Eq.(14) in (13) and equate coefficients of ^ yields 



(14) 



&i = -!-(a±2>/=£) (15) 
la 

When it comes to the discussion of the behaviour of x at infinity, one expands x as: 



A Ai 

X = a + - + — (16) 

y y 2 

and find A as 



V 2 

X = ± T J 1W (17) 
One can see that the behaviour of x is for large y. Hence, 



b 1 +n = X (18) 

In order to find the wavefunction, x(y) can be written as the sum of the Laurent expansions around 
different singular points, plus a constant C\. Hence 



where P n (y) is a n th degree polynomial. Substitute Eq. (19) in (13) and get 



P n P n \y J V 2/ 2/ a aWVi/ 

For large y one can find Ci = The wave function in terms of x can be written by using eq. (19) 
and (7) as 



^) =e *H/|Vf (21) 
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In Eq. (21) the correct value of C\ is used as C\ = — ^ because of the condition for the wavefunction 
which is known as y — > oo, ip(y) — > 0. It is seen from Eq. (15) that bi has two values and no particular 
value has been chosen. Using Eq. (18), the energy eigenvalues for any n-th state become, 



E n — 



a 



-(2ra + l) + 



V, 



If we use Eq.(22),(15) and C\ = —1/a in Eq.(20), it becomes 



(22) 



vKiv) + 



+ l-y)P^y)+nP n (y) = 



'2a 

which is a Laguerre differential equation. Therefore, the wavefunction is obtained as 



(23) 



MV) = C n e-i y ( ~ n± ^ ] L 2 n a ^(y). (24) 

where C n is a normalization constant and Ln"^ (y) are Laguerre functions. The wavefunction satis- 
fies the boundary condition that is y — > oo, ip(y) — > 0. 

3.1 Non-PT symmetric and non-Hermitian Morse Potential 

In equation (10), if the potential parameters are defined as V\ — {A + iB) 2 , V 2 = (2C + 1)(A + iB) 
and a = 1, then the potential becomes [19], 

V(x) = {A + iB) 2 e' 2x - (2C + 1)(A + iB)e~ x (25) 
where A, B and C are arbitrary real parameters and i = >/— T. The QHJ equation is 



p 2 - ip' - [£ - (A + iB) 2 e~ 2x + (2C + 1)(A + i^e - *] = (26) 

Using the transformation in the form of y = (A + iB)e~ x in the Eq. (26), then using = iy<p and 
X = + ^, Eq. (26) becomes 

7 







4y2 y. 

As it is seen from Eq. (27), x has a pole only at y = and for y = define % for the Eq.(27) as, 



(27) 



6i 

X = — + a + a x y 

y 



Using the Eq. (28) in (27), b\ is found as 



Following the same procedure that is given in section 3, energy is found as 



E n = -{n - Cf. 



(28) 



(29) 



(30) 



Now choose the potential parameters in Eq.(10) as V\ is real and V 2 = A + iB, the Morse potential 
can be written in the following form 



V(x) = V ie - 2iax - (A + iB)e 



(31) 



and following the same procedure, the energy is obtained as 



E n = a 4 



(n + 1/2) 



A + iB 



2aJ\ - Vx 

According to Eq. (32), the spectrum is real in the case of Im(V 2 ) = 0. 



(32) 



3.2 PT symmetric and non-Hermitian Morse Potential 

When a = ia and V±, V 2 are real, the Morse potential becomes 
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V(x) = V ie - 2iax - V 2 e 



(33) 



Thus, the energy eigenvalues are obtained as 



E n = a 4 



Vo 



(34) 



2 2a y J\-V 1 \ 

If we take the parameters of Eq.(39) as V\ = —oo 2 , V 2 = D and a = 2 then, corresponding eigenvalues 
for any n-th state are obtained as 



E n = (2n + 1 + £-) 2 



are consistent with the results [10-11,19]. 



4 Poschl-Teller Potential 



The general form of the Poschl-Teller potential is 



V(x) = -W 



-2ax 



(1 + qe~ 2ax y 



The QHJ equation is given as 



(35) 



(36) 



— 2ax 



p 2 -ip - \ E + AV Q - 



If we take y = ±.i^fqe~ ax and use the transformations as p = 
equation turns into 



(37) 



—iaycj) and x — 4> + ^> the QHJ 



X+X + 



1 1 

Ay 2 a 2 y 2 



E 



AVr 



o y 



q {l-y 2 ) 2 



= 0. 



(38) 
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As it is seen from Eq.(38), x has poles at y = and ±1. % is expanded at y = and bi is found 



2a 

At y = 1, one can expand % as 



h = — (a±2 v /z £) (39) 



6' 

X= -^— + a' + a 1 (l-y) (40) 

i - y 



and b[ is found as 



&i = ^-(?«± V«V + 8gF ) (41) 
At |/ = — 1, one can expand x as 

X= Y^ + «o + «i'(l + y) (42) 

Substitute Eq. (42) in (38), to obtain 6^ = 6^'. One can look at the behavior of x at infinity with 
expanding x as 



A A ^ , . 

X = A + - + — (43) 

v y 

From Eq. (43) and (38), A is found as 



A = — (a±2^E) (44) 
2a 

and behavior of x is 61+6 i+ fc i +2ra f or large y. Hence, 

\ — bi + b' 1 + b\ + 2n (45) 
In order to find the wavefunctions, x is written as 
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x= h + ^ + j^ + m +c , (46) 

y i-y l + s/ inU/j 

Substituting Eq. (46) in (38), C 2 can be found as C 2 = for large y. The wavefunction can be written 

as 



(/(^+tV4+BM 

If we look at Eq. (47), fti and ^ have two values and both of them gives appropriate results for 
energy spectrum and wavefunction, no particular value has been chosen. Thus, residues are given as 
bi = ^(a ± 1\J — E) and 6^ = b" x = ^(qa ± \fd 2 q 2 + 8gV ). Finally the energy is obtained by 
using Eq. (45) as 



a 2 /\ , x 8V? 



E n = -— (2n+l)±Wl + — £ (48) 



4 \ v V Q ,Q;2 

Using Eqs. (47,48) and (38), one can find the wave function as 



My) = n y -(-^)±7 (1 _ ^i^) ^-i (y) (49) 

where N is a normalization constant, 7 = ,/l + P n " 2 2 '" 2 2 (y) stands for Jacobi polynomials 

V 

and vi = If we look at Eq.(49), there are three cases for physical solutions because of the 

wavefunction that satisfies the boundary condition as y — > 00, ip(y) — > 0. If — (n — 1/2) ± 7 < 
andl±7 > 0, it should be 1 ± 7 > | - (n - 1/2) ± 7]. If -(ra-l/2)±7< 0andl±7 < 0, there 
is no restriction for the parameters and there are physical solutions in this case. The last case can be 
defined as; if — (n — 1/2) ± 7 > 0, 1 ± 7 > | — (n — 1/2) ± 7] for appropriate solutions. 



4.1 Non-PT symmetric and non-Hermitian Poschl-Teller cases 

Here, V and q are complex parameters V = V 0R + iV 0I and q = q R + iqj but a is a real parameter. 
Although the potential is complex and the corresponding Hamiltonian is non-Hermitian and also non- 
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PT symmetric, there may be real spectra if and only if V 0I q R = V 0R qj. When both parameters V and 
q are taken pure imaginary, the potential turns out to be, 



V(x) = -Wo 



2 qe - 4ax + i(l - q 2 e - 4ax ) 



(50) 



(1 + g 2 e- 4 ^) 2 

For simplicity, we use the notation V and q instead of V 0I and qi. In this case, we get the same energy 
eigenvalues as in Eq.(48). If q is an arbitrary real parameter and V =>- iV also a =>- ia completely 
imaginary, the potential becomes 



V{x) = -W c 



(1 — q 2 )sin2ax + i(2q + (1 + q 2 )cos2ax) 



[1 + q 2 ) 2 + Aqcos2ax(l + qcos2ax + q 2 ) 
and the corresponding energy eigenvalues become 



(51) 



E 



a 

T 



2n + l + \ 1 + 



8V 



(52) 



For simplicity, let us take all three parameters a, q, V purely imaginary. Then the potential takes the 
form 



V(x) = -4V 



and the energy becomes 



(1 + q 2 )sin2ax + 2q + — q 2 )cos2ax) 
(1 + q 2 ) 2 + 4g 2 (l - cos 2 2ax) + 4g(l + q 2 )sin2ax 



(53) 



a 



2n + 1 + ^-Ja 2 q 2 + (l + q 2 )V 
2aq v 



(54) 



4.2 PT symmetric and non-Hermitian Poschl-Teller cases 

We choose parameters V and q real and also a = ia. Then, the potential turns into 
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V(x) = -W 



(1 + q 2 )cos2ax + 2q + i(q 2 — l)sin2ax 



(55) 



(1 + q 2 ) 2 + 4qcos2ax((l + qcos2ax + q 2 ) 



and corresponding energy spectrum is 




2 



(56) 



5 Conclusions 



We have applied the PT-symmetric formulation to solve the Quantum Hamilton- Jacobi equation for 
Morse and Poschl-Teller potentials in both real and complex forms. We have obtained the energy 
eigenvalues and the corresponding wave functions for different forms of these potentials within Quan- 
tum Hamilton- Jacobi formalism. The real energy spectra of the PT-/non-PT- symmetric complex val- 
ued non-Hermitian potentials have been obtained in case the potential parameters are restricted. It is 
also shown that the QHJ formalism is a good approach to obtain eigenfunctions and energy eigenval- 
ues for a class of exponential type potentials discussed here within PT symmetric frame. As a result, 
we have pointed out that our exact results of complexified general Morse and Poschl-Teller potentials 
may increase the number of applications of complex Hamiltonians with real energies in the extensive 
study of different quantum systems within the flexible Quantum Hamilton-Jacobi approach. Finally 
we should state that this work is the first application on the study of PT-symmetry for the Quantum 
Hamilton-Jacobi approach. 
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